PHYSICAL REVIEW LETTERS Submitted 
Dissipative Scaling Functions in Navier-Stokes Turbulence: Experimental Tests 
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A recent theoretical development in the understanding of the small-scale structure of Navier- 
Stokes turbulence has been the proposition that the scales r]n{R) that separate inertial from viscous 
behavior of many-point correlation functions depend on the order n and on the typical separations 
R of points in the correlation. This is a proposal of fundamental significance in itself but it also 
has implications for the inertial range scaling behaviour of various correlation functions. This 
dependence has never been observed directly in laboratory experiments. In order to observe it, 
turbulence data which both display a well-developed scaling range with clean scaling behaviour and 
are well-resolved in the small scales to well within the viscous range is required. We have analysed 
the data of the experiments performed in the laboratory of P. Tabeling of Navier-Stokes turbulence 
in a helium cell with counter-rotating disks, and find that this data satisfies these criteria. We have 
been able to find strong evidence for the existence of the predicted scaling of the viscous scale. 



' The famihar approach to the statistical theory of 
, Navier-Stokes turbulence concentrates on the proper- 
ties of two-point differences of the Eulerian velocity field 
, u(x, t) and their moments, termed structure functions: 



S^{R) ^ {\u{r + R) - u{r)Y') 



(1) 



In isotropic homogeneous turbulence, these structure 
functions are observed to behave as a power-law in R, 
Sn {R) ~ , with scaling exponents that may be 
universal. This scaling holds within a range of scales be- 
tween the outer scale L determined by the system size 
or the forcing, and some inner scale 77 determined by 
the viscosity, below which the velocity field is essentially 
smooth. In this regime 5„(i?) ^ i?". The usual defini- 
tion of the viscous scale was established by Kolmogorov 
from the balance of the viscosity v and the mean energy 
flux e, according to 77 ~ iRe^^^'-^^''^'' . The observation of 
intermittency in turbulence is contrary to the relatively 
straightforward picture which follows from Kolmogorov's 
assumption of a single length scale, and suggests that a 
more complex situation applies. We begin by a review of 
the derivation of the functional behaviour of the viscous 
scale (found in more detail in |2|). 

One possible approach to the question of cross-over 
scales follows a recent theoretical trend to concentrate 
on the more general simultaneous many-point correla- 
tion functions J-"„ of velocity differences rather than on 
two-point quantities only. These are defined in terms of 
two-point differences w{r,r',t) = u{r',t) ~ u{r,t), as 



{w{ri,r[)w{r2,r2) . ..w{rn,r'J) , 



(2) 



where (•) denotes averaging, and all coordinates are dis- 
tinct. Time labels have been dropped as only simultane- 
ous correlations will be considered. Homogeneous scaling 
here means that 

^„(Ari,Ar'i,...,Ar;) = A^":F„(ri,r;,...,0, (3) 



with C„ the scaling exponent. Taking the time deriva- 
tive of using the Navier-Stokes equations to evaluate 
each du/dt and considering the stationary state where 
dlF/dt = 0, one derives the following statistical balance 
equation: 



Vn{ri,r[; ...;r„,r;j = vj„{ri,r[; ...;r„,0 



(4) 



Here the term arises from the nonlinear interaction 
term and may be written as 
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X (u;„i(ri,r'i) . . . L'^{rj,r'.,r) . . . u'a„ (r„, r^,)} , (5) 



L^(r„r',r) 
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^w^ivj - r,rfc)V 
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In the above Paj/iir) is the projection operator. The 
RHS with coefficient the kinematic viscosity, results 
from the viscous term and is defined 



3 = i 



X {wa,{ri,r\) . ..WaArj,r'A . . .7«a„(r 



(6) 



This equation provides the means to determine the scale 
of the viscous range. The balance equation expresses the 
competition between the small-scale viscous effects and 
the interesting non- linear dynamics, and intuitively, the 
viscous scale should be the scale at which the two effects 
become comparable. This raises a rather subtle question. 
As there is a balance of the two terms, it appears that the 
scaling properties of the correlators must be determined 
by the viscous term. However one believes that the prop- 
erties of the inertial range quantities are independent of 
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the details of the viscous range. This apparent paradox 
can be understood if one considers the separations of the 
coordinates of the correlation functions. If all separa- 
tions are in the inertial range, there are no small-scale 
quantities, and the contribution of Jn will be negligible, 
leaving a homogeneous equation I?„ = 0. Non-trivial 
scaling can arise from special solutions for the terms in 
the sum in P„ that exactly cancel one another. Now as 
some coordinates in the correlation functions approach 
one another, the gradients in Jn will begin to show their 
effect: these pick up the smallest separation rmin, intro- 
ducing a factor of l/r^^j^. As Tmin 0, this term is 
no longer negligible; the scaling solution of the homoge- 
neous inertial range equation will no longer be valid and 
one obtains the smooth viscous result. 

Thus one wishes to estimate the terms of the balance 
equation both in the inertial range and in the limit where 
some coordinates approach one another, in order to ob- 
serve this crossover and estimate its scalelength. 

It can be shown that in the case where all separa- 
tions are of order R one may evaluate V simply as 



{R)IR 



(7) 



This can be demonstrated by proving that the integral in 
(||) converges in both limits, so that the typical evaluation 
at R is correct. As points in the correlation approach one 
another, this evaluation remains valid; but cancellations 
between terms will no longer occur. 

The second term Jn can be estimated directly as 
Jn{R) ^ Sn{R)/R^ ■ In the limit when some separa- 
tion Ty — > 0, this evaluation is replaced by Jni^ij', R) ~ 
^n{i"ij ) R) / fij 1 where Tn{rij]R) is shorthand notation 
for Tn with an overall typical separation R and some 



pair of coalescing points of smaller separation 
the balance equation gives 



rlSn+i{R)lvR 



Now 



(8) 



This gives the evaluation of Tn {r-ij ; R) for a small sepa- 
ration in the viscous regime. 

Now we wish to compare this with an evaluation for 
[fij ; R) when the small distance is still in the inertial 
range. To do so we invoke the fusion rules derived in . 
These rules predict the behaviour of multipoint correla- 
tion functions as some pairs of coordinates approach one 
another, or "fuse". The essential result concerns a cor- 
relation of n pairs of points JF„ where p pairs of coordi- 
nates ri,r[ . . .rp,r'p, (p < n) of p velocity differences co- 
alesce, with typical separations between the coordinates 
~ ''il for i < p, and where all other separations 
are of the order of i?, r <^ R <^ L. Let us denote such a 
correlation as J-'^n'^ ■ In a homogeneous isotropic scaling 
system, the fusion rules predict 



(9) 



•^pi'^l: ^1 1 ■ ■ ■ I '^pi '^p)^n,p {fp+li i ■ ■ ■ I '^ni ^n) 



where j^p is a tensor of rank p associated with the first p 
tensor indices of and it has a homogeneity exponent 
(p. The (n-p)-rank tensor ^n,p{rp+i,r'p_^i; ■ ■ . ;r„,a;^) 
is a homogeneous function with a scaling exponent ^„ — 
Cp, and is associated with the remaining n — p indices of 
JF„. In terms of the scaling of structure functions, this 
can be expressed for p points coalescing to a distance r 
and all other points with typical separation R as (abbre- 
viating the coordinate dependence of ^n ^) 



TiP\r;R)^Sp{r)Sn{R)/Sp{R) 



(10) 



In the special case that p = 1, due to the vanishing of 
the average of a single difference in isotropic turbulence, 
the leading order result is 



Ti'\r-R) ^ S2{r)SniR)/S2{R). 



(11) 



Applying this result to the correlation we previously de- 
noted J^ni^ij] R) one obtains 



TniUj-.R) - S2{nj)Sn{R)/S2{R). 



(12) 



Thus we have an inertial range and a viscous range 
evaluation of J-n(rij; R). Let us take the scale r]n to be 
that at which the two evaluations coincide. Balancing in 
the two-point case one recovers the Kolmogorov estimate, 
772 ^ LRe"^^*-^"^^-*. For other values of n one finds 



Vn{R) = m 



R 



Cn + Cs - C 



n+1 



C2 



C2 



(13) 



In order to test this proposition experimentally, we will 
consider direct measurements of the function Jn{R)- To 
make a comparison with the one-dimensional data ob- 
tained from experiments we take a form defined by 



J„(p;i?) = (vXr) [w{r,r + R)Y 



■ R/R. (14) 



For discrete data, the Laplacian operator in (|lj) is 
taken to be a second order finite difference of longitudinal 
components of the velocity, 

V^u(r) = [w{r, r + p) - w{r, r - p)] • p/p^. (15) 

From the discussion above, one expects a different scal- 
ing for p above and below the dissipative scale. For p in 
the inertial range, the estimation of (^ is applicable and 
one predicts 

Jn{p; R) = CnJ2{p)Sn{R)/2S2{R) , p » 77 . (16) 

where C„ is an i?-independent dimensionless constant 
which may have n-dependence. However, for p in the 
viscous regime, 

Jnip; R) = CnJ2{p)Sn+l{R)/S3{R) , P « 7? , (17) 
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FIG. 1. Log-log plot of the structure functions SniR) as a 
function of R for n = 1 — 10. 
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FIG. 2. Log-log plot of the normalised function 
J4,{R) / J2{R) as a function of R for p =1,3,7,11 and 29, rep- 
resented by □, +, X, * and o respectively. The data are 
compensated in the upper plots by the inertial range fusion 
rule prediction S4,{R) / S2{R) and in the lower by the result in 
the viscous regime S5{R)/Sz{R). 

where (?„ is some other coefficient. One can show that 
J2 is equal to the mean dissipation (|Vu(a:)p), and is 
thus expected to be i?- independent. The explicit prefac- 
tor containing p is included in J2; we will consider only 
the R dependence resulting from the scaling of ([l^). 

These predictions are tested in data obtained by F. 
Belin and H. Williame in the laboratory of P. Tabeling 
at Ecole Normale Superieure; see for example (|,^. The 
data are time signals of the velocity field taken from a 
low-temperature cell of helium gas enclosed in a cylin- 
der and driven by counter-rotating disks. The helium is 
maintained at a constant temperature around 5K and at 
a controlled pressure. The recordings were made using a 
hot-wire probe consisting of a 7/xm carbon fibre coated 



with evaporated gold apart from an active area of size 
of order lO^m. The frequency response of the probe can 
range between 10 and 50 kHz. The data had very long ac- 
quisition times, containing up to 30 million samples. The 
resulting statistics are well-resolved and stationary. The 
Taylor microscale Reynolds number and the Kolmogorov 
microscale 77 were determined through the usual proce- 
dure of surrogating time for space (by Taylor's hypothe- 
sis), and data is available for a range both of R\ and of 
minimum resolved lengthscale r/r]. We selected data sets 
according to the small-scale resolution, and although the 
R\ was not extremely high, a distinct inertial range is ev- 
ident. The data presented here has a minimum resolved 
distance r/rj of 1.18 and R\ of 418. 

In Fig. we present the structure functions Sn{R) as 
a function of i?. In all figures, spatial separations have 
units of sampling times, and the velocity is normalised 
by the RMS velocity. This figure shows that we have one 
and a half decades of "inertial range" (between, say 10 
and 500 sampling units) and that the lengthscales below 
10 units are smooth and well-resolved. The initial loga- 
rithmic slope of the 71th structure function at 1 unit is 
close to n (deviating successively more for higher order n, 
as would be expected). There is reasonably well-defined 
scaling behaviour up to order 10. 

Our aim is to try to expose the postulated cross-over 
in scaling behaviour in R of Jn {p', R) as a function of p. 
We have calculated the correlation functions Jn(p; R) for 
several values of p from the minimum distance of 1 unit 
up to a value well into the inertial range. Note that the 
difference in scaling that is expected is rather small; one 
expects the scaling exponent of J„(p;i?) to cross over 
from (n — C2 to C,n+i — Cs) which for the usual values 
of scaling exponents obtained in turbulent experiments 
gives a difference of the order of 0.15 for n = 4 to 0.2 for 
n = 8. Thus we do not present the results in terms of 
calculated exponents, as one cannot justifiably separate 
values of this closeness on the basis of exponents calcu- 
lated on a limited inertial range. Instead we will examine 
the function as a whole. 

In Figs. 2-4 we display the results. The three figures 
show Jn{p', R) for a single value of n, for n = 4, 6 and 8. 
For each n there are results for five values of p, p =1, 
3, 7, 11 and 29. The figures each show two sets of 
data, one in which the calculated J„s have been compen- 
sated by the inertial range prediction J2{p)Sn{R) / S2{R) 
(the upper set of functions), and the second showing the 
same data compensated by the dissipative range result 
J2ip)Sn+iiR)/ S3{R). Hence in the upper set we expect 
to see that for inertial range values of p, the resulting 
plots are constant in R in the inertial range. In principle 
as there is no knowledge of the coefficient C„, the value 
of the constant C„ can be different for different n. (It is 
trivially 1 for n = 2.) We hope to see that the dissipative 
range scaling is a better fit as p ^ 0. 
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FIG. 3. Log-log plot of the normalised function 
Je (/?)/ J2 (i?) as a function of R for p =1,3,7,11 and 29, repre- 
sented by □, -f , X, * and o respectively, compensated in the 
upper plots by Se{R)/S2{R) and in the lower by St{R)/S3{R). 



ing show a distinct indication that there is a tendency 
toward this slope. The effect becomes more apparent for 
larger n as the separations between the two scalings be- 
comes larger. These functions of course also show larger 
statistical fluctuations. 

One should note that the plots for p = 1 and p = 3 
are almost identical. This shows that in fact the data is 
only resolved well to p = 3, and no further information is 
gained in the subsequent refinements of scale. Thus the 
minimum lengthscale that can be considered as resolved 
is in fact of order It is therefore not surprising that a 
clean scaling of (^7|) is not precisely observed. Nonethe- 
less the trend in that direction is clearly visible. 

We have been able to find convincing experimental ev- 
idence that the viscous scale of n-point multipoint cor- 
relation functions is an anomalous scaling function. We 
have verified the inertial range fusion rules and given ev- 
idence that the small scale structure behaves according 
to the theoretical predictions. 
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FIG. 4. Log-log plot of the normalised function 
Js{R) / J2{R) as a function of R for p =1,3,7,11 and 29, repre- 
sented by □, -f , X, * and o respectively, compensated in the 
upper plots by Ss{R) / S2{R) and in the lower by Sg{R) /Sz{R). 

The upper sets of plots in each figure show that for 
the inertial range values of p, the inertial scaling (^6|) is 
very well realised. This scaling has been previously ob- 
served in turbulence data 0| but the agreement in this 
data is more impressive: it has smaller fluctuations, and 
the agreement continues into the viscous scales, which 
has not previously been seen to be the case. Comparing 
between the figures, for different values of n one finds 
that all coefficients C„ are very near to 1. Comparing 
different values of p, there is a continuous dependence on 
p in the functional behaviour of J„ (p; R) . There is a clear 
deviation from the inertial range scaling as p decreases, 
and the smallest value of p shows a small but definite 
slope. The plots corrected by the dissipative range scal- 
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